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We have analyzed the consequences of introducing the local density approximation combined with an effective 
nuclear momentum-dependent potential into the CC quasi-elastic neutrino-nucleus scattering. We note that the 
distribution of the recoil nucleons momenta becomes smooth for low momentum values and the sharp threshold 
is removed. Our results may be relevant for Sci-Fi detector analysis of K2K experiments. The total amount of 
observed recoil protons is reduced because some of them remain bound inside the nucleus. We compare theoretical 
predictions for a probability of such events with the results given by NUX+FLUKA MC simulations. 



1. INTRODUCTION 

In recent years there has been growing interest 
in the studies of neutrino interactions at energies 
of a few GeV 1 . It was motivated by the need 
for more precise measurements of neutrino oscil- 
lation parameters (#13 in particular). This entails 
deriving the best description of interactions with 
free nucleons, and then incorporating nuclear ef- 
fects. From the point of view of the Monte Carlo 
codes most (or all) nuclear effects are described 
with numerical packages |2] but it is enlightening 
how many of these effects can be presented in an 
analytical form. 

We investigated the process n — ► M P 
with the target neutron bound inside the nucleus. 
Computations of nuclear effects were based on the 
Fermi gas model which is known to work well 
in the above mentioned energy region But 
this approach is not completely satisfactory. For 
example its simple form leads to the conclusion 
that the ejected nucleons can only have momenta 
higher than the chosen value of the Fermi momen- 
tum kp However there seems to be no phys- 
ical reason why lower values of momenta should 
be forbidden. 

A possible solution to this problem is to intro- 
duce a local density approximation . In fact the 
density of nuclear matter is not uniform and ac- 
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cordingly it is possible to introduce the concept 
of local Fermi momentum kp(r). Since the in- 
teraction can take place in a region where Fermi 
momentum is arbitrarily low, the distribution of 
momenta of recoil nucleons becomes smoother. 

Another solution can be based on a different 
effect. The target nucleons are not free but they 
are bound inside the nucleus. The binding en- 
ergy is smaller then the typical values of energy 
transfer but it may cause interesting effects. Mo- 
mentum dependent optical potential can be ob- 
tained even from the simplest versions of nuclear 
mean field theory 0. In the covariant approach 
based on the Dirac equation one has to distin- 
guish contributions from scalar and vector com- 
ponents of nucleon self-energy in nuclear matter. 
In this framework the electron-nuclei scattering 
was discussed in [7] with mean fields taken from 
More recent results on the self-energy can 
be found e.g. in jJJ] where derivations used the 
G-matrix approach. Our investigation is based 
on other computations This choice was dic- 
tated by two reasons. Firstly we were able to 
reproduce an explicit formula for the potential 
as a function of density (local Fermi momentum) 
and nucleon momentum. Secondly we wanted to 
enable comparisons with the numerical results of 
other authors who used the same potential 

We wanted to describe both local density and 
potential effects and for that we needed an ana- 
lytical form of the potential. So we have derived 
a simple analytical form of the potential depen- 
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dent on two variables: the Fermi momentum and 
the nucleon momentum. We incorporated this 
potential into a Monte Carlo generator of events 
and obtained a spectrum of the ejected nucleon 
momentum. The simulations were performed for 
three target nuclei: oxygen, iron and argon; these 
are possible targets in the neutrino experiments. 
The results in all three cases were very similar be- 
cause the shapes and the density profiles did not 
differ significantly. 

We investigated also the possibility that the ef- 
fective potential leads to reduction of the number 
of the produced protons that escape from the nu- 
cleus. We compute the fraction of the total cross 
section which can be interpreted as correspond- 
ing to an excited nucleus in the final state. The 
same effect is predicted by other MC generators 
e.g. by NUX+FLUKA. 

We have also looked for the effects related 
to the proton-neutron asymmetry inside the nu- 
cleus. It turned out that the introduction of sep- 
arate values of the Fermi momentum for proton 
and neutron gases left all the plots virtually un- 
changed. 

Our results can be useful in near detector anal- 
ysis of neutrino interaction events in K2K where 
Sci-Fi detector registers tracks of ejected protons. 
The reconstruction procedure requires the parti- 
cle to cross at least 3 planes of scintillating fiber 
thus the protons would need a momentum greater 
than 500 MeV. Therefore in the data analysis it is 
important to have MC which correctly describes 
the shape of the distribution of recoil protons mo- 
menta. 

We performed MC simulations assuming the 
neutrino energy profile expected at Sci-Fi detec- 
tor ■ We investigated to what extent the spec- 
trum of the recoil proton momentum is influenced 
by the introduction of the LDA and the effective 
potential. 

We restricted ourselves to a simplified dynamic 
model without the effects caused by the RPA cor- 
relations. We also excluded form our analysis the 
effects due to the final state interactions in nuclei. 

The spectrum of the ejected nucleon momen- 
tum was studied earlier in the context of NC re- 
actions where it is the only observable quantity. 
The outgoing nucleons with higher momentum 



(above the Chcrcnkov threshold for SK i.e. big- 
ger than about 1.07 GeV) were studied recently 
in |13| . In an earlier study |14j an average binding 
energy was used to describe kinematics of the pro- 
cess and interesting predictions for angular distri- 
butions of recoil nucleons were proposed. Pauli 
blocking was not imposed and arbitrarily low val- 
ues of the momenta of ejected nucleons were ob- 
tained. Another approach [TSj used Pauli block- 
ing and forbode recoil nucleons of kinetic energy 
lower than about 27 MeV. It seems clear that any 
realistic model of interaction should incorporate 
Pauli blocking and consequently also a mecha- 
nism to remove the nonphysical threshold in the 
ejected nucleon momentum spectrum. 

2. MODEL 

Our Monte Carlo generator simulates quasi- 
elastic neutrino interactions with basic dynamics 
introduced according to |lt>| . 

The events are obtained in the following man- 
ner: 

• The neutrino energy E v is chosen as either a 
fixed value or generated according to some 
beam energy profile. 

• The Fermi momentum is established using 
global or local Fermi momentum scheme: 

— In the global scheme the Fermi mo- 
mentum is fixed. 

— In the local scheme the region in the 
nucleon where the interaction is going 
to take place is selected first. Then the 
Fermi momentum is calculated based 
on the nuclear density in this region. 

The actual target momentum is chosen at 
random from the Fermi ball of that radius. 

• The nucleon neutrino pair is boosted to its 
center of mass frame (CMS) where the di- 
rection of the scattering is taken to be ran- 
dom but it will be weighted later on. 

• While keeping this direction fixed in the 
CMS frame various values of the outgo- 
ing nucleon momentum are tested and by 
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means of the bisection algorithm the value 
which brings about the energy conservation 
is chosen. The energy is always evaluated 
in the LAB frame and takes in the account 
the momentum dependent potential of the 
nucleon. 

• The outgoing nucleon is Pauli blocked if its 
momentum in the lab is smaller than the 
local fermi momentum. 

• The exit of the nucleon from the nuclear 
matter is simulated by diminishing the val- 
ues of its momentum from the above calcu- 
lated value pn (inside nucleus) to the value 
p' (outside nucleus) calculated from the en- 
ergy conservation condition 



V{p N )+E k {p N ) = E k {p') 



(1) 



If there is no solution then the proton is 
unable to leave the nucleus and the nucleus 
becomes excited. The justification for this 
equation comes from the fact that there is 
no nuclear potential outside the nucleus. 

• The value of cross section is calculated ac- 
cording to neutrino energy in the target rest 
frame. A correction of little significance 
(less than 1 % effect on the cross section) 
due to nonlinearity of the neutrino and tar- 
get nucleon momenta is taken into account. 
The weight of the event is proportional to 
the nuclear density, differential cross section 
(with correction), the boost jacobian, the 
bisection algorithm jacobian. 

The target is treated as a collection of nucle- 
ons distributed in space according to the density 
profile determined by electron-nucleus scattering 
experiments. For & O w we adopted the harmonic 
oscillator model in which the nuclear density is 
given by 

p 16 °(r)= Po exp(-r 2 /i? 2 )(l + C^) (2) 

where R = 1.883 fm, p = 0.141 fur 3 , C = 1.544, 
and po is a normalization constant defined by the 
condition: 




Figure 1. The local Fermi momentum r depen- 
dence for 8 16 , ls Ar 40 and 2 6^e 56 . 



For isAr 40 and 2oFe 56 we use the two param- 
eter Fermi model and write the density as: 
Po 



Ar,Fe 



(r) 



r-C 
Ci 



1 + exp 

with the following parameters: 



(4) 







C[fm] 


Ci[fm) 




0.176 
0.163 


3.530 
4.111 


0.541 
0.558 



The local Fermi momentum is determined by 
the density profile according to: 



k F (r) 



3 3n 2 p(r) 



(5) 



d 3 rp °{r)=A. 



(3) 



In fig. ^ we show the local Fermi momentum 
dependence on r for Oxygen, Argon and Iron. In 
the case of nonsymmetric nuclei one can also in- 
troduce separate local Fermi momenta for protons 
and neutrons: 

Fin 

where A, Z, and N are the atomic number, the 
number of protons, and the number of neutrons 
in the nucleus respectively. 
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Figure 2. Momentum dependent potential 
V(kp,p) for 3 values of Fermi momentum (see 
formula (7)) compared with three original plots 
taken from 



The average value of fcjr is calculated as: 

/ k F (r) r 2 p nucleus (r) dr 



/ ^nucleus ^ 



j.2 pnucleus ^y^j 



(7) 



We get the following values: 

< fc° >= 199 MeV, 

< kj r >= 217 MeV, 

< fc£ e >= 217 MeV. 

Using tresults from ^U] we find an analytic 
form of the real part of the optical potential as 
(see fig-HJ: 



V(k F , P ) 



(ak F ) 2 (k F + b) 
d 3 k f + e 3 p 2 /k f + p 4 ' 



(8) 



where kF, P, and V(kp,p) are given in MeV. The 
fitted values of the parameters are: 
a = 206 MeV, b = 582 MeV, c = -322 MeV, 
d = 422 MeV, and e = 289 MeV. 

The formula © fulfills the following criteria: 

• V is negative 



• V reproduces essential features of plots from 

ma, 

• V is monotonously increasing, 

• For higher values of momentum, V quickly 
approaches zero, 

• For small values of momentum, V is pro- 
portional to p 2 as in the limit as p — > on 
general ground one expects that 



V(p) ~ V + 



P 



2M* 



(9) 



There is some amount of uncertainty in the 
reconstructed potential which is inherited from 
the original computations in which contributions 
from only finite number of harmonics were in- 
cluded. 
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Figure 3. Total cross section for quasi-elastic 
scattering on iron nucleus per nucleon. 



In numerical computations we compare three 
cases: 

(a) Fermi gas with global kp = (k F ucleus ^ 

(b) Fermi gas in local density approximation 

(LDA), 



5 





0.11 r 




0.1 - 




0.09 - 




0.08 - 


c 


0.07 - 


o 




o 


0.06 - 


ffl 


LL 


0.05 - 




0.04 - 




0.03 - 




0.02 - 




0.01 - 



Potential in this paper 
Nux+Fluka Monte Carlo 



0.5 1 1.5 2 

neutrino energy (GeV) 



2.5 



> 
o 
CD 



12 



10 







Kp = 199 MeV — 






"\ LDA 






bQA and potential - . 









0.2 0.4 0.6 0.8 1 

recoil nucleon momentum (GeV) 



1.2 



Figure 4. A fraction of the total cross section 
corresponding to events in which proton remains 
bounded in excited nucleus. 



Figure 5. Quasielastic neutrino-oxygen scatter- 
ing: recoil nucleons momentum distribution at 
E v = 700 MeV. 



(c) Fermi gas with in the local density approxi- 
mation with momentum dependent nucleon 
potential ©. 

3. RESULTS 

The total cross sections for quasi-elastic inter- 
action for all three cases are plotted in fig. [3J 
While changes introduced by the LDA are mi- 
nor, the effective potential reduces the total cross 
section by a few percent. Since the produced pro- 
ton is affected by momentum depended potential 
it may not have enough energy to leave the nu- 
cleus. In such a case in the final state the nucleus 
is excited and there is no ejected nucleon. This 
happens always if 

V( PN ) + E k {p N ) <M. (10) 

The probability of such a process is shown in 
n g-El F° r comparison we calculated the probabil- 
ity of analogous events in the NUX+FLUKA MC 
generator. It is very interesting that the shapes of 
the curves are in both cases almost identical. The 
probability for a bound nucleon in the final state 
decreases with neutrino energy in both genera- 
tors and at energy of about 1 GeV becomes flat. 



However, the effect predicted by our Monte Carlo 
generator is much bigger (about three times). 

The predicted distributions of the ejected nu- 
cleon momentum are shown in fig. [5] and |5| The 
plots are normalized as differential cross sections 
da I dp' . In all cases the neutrino energy is fixed 
at E v — 700 MeV. The changes observed due 
to the LDA and the effective potential are very 
similar for all three analyzed targets so we do not 
include a plot for iron. 

In the case (a) there are obvious sharp thresh- 
olds at (k F ucleus ). 

In the case (b) the distributions become 
smoother. They differ only in the region of nu- 
cleon momentum below ~ 300 MeV; whereas, for 
higher values of nuclear momentum, LDA intro- 
duces no significant differences. 

Finally, in the case (c) there is higher probabil- 
ity for the outgoing nucleon to have momentum 
lower than kp. The values close to zero become 
accessible as well. The differential cross section 
rises almost linearly for nucleon momenta smaller 
than 450 MeV then it slowly bends reaching max- 
imum at about 600 MeV. For momenta above 
800 MeV the plot is only slightly changed com- 
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Figure 6. Quasielastic neutrino- argon scattering: 
recoil nucleons momentum distribution at E v = 
700 MeV. 



Figure 7. Effect of different Fermi momenta for 
proton and neutron Fermi gases in quasielastic 
neutrino-iron scattering for ejected momentum 
distribution at E v = 700MeV. 



pared to the cases (a) and (b). 

The fig. illustrates the relevance of the intro- 
duction of separate values of Fermi momentum 
for neutron and proton Fermi gases. The Fermi 
momentum for neutrons is larger and since the in- 
teraction takes place on neutrons, it is clear that 
the cross section should be slightly larger. It is 
indeed the case but the effect is very small and 
the difference of the plots is hardly noticeable. 

For completeness in fig.[S]we present how both 
discussed effects influence the energy transfer 
spectrum for the iron nucleus (plots for oxygen 
and argon are very similar). There is only a 
small difference at the peak between the differ- 
ential cross section plots for the cases (a) and (b) 
the plot for the LDA being slightly lower. In the 
case (c) we observe a significant change in the 
shape of the plot. For lower values of the energy 
transfer the differential cross section is substan- 
tially reduced and for higher energy transfers it 
is slightly increased. The allowed kinematical re- 
gion is also modified. 

In fig. we show the ejected proton momentum 
distribution obtained for neutrino beam with the 
energy profile identical to that predicted for the 



K2K near detector. We note that the LDA and 
the potential effects remove the sharp threshold 
and give a smooth plot at low values of the mo- 
mentum. The number of nucleons with momen- 
tum between and 700 MeV is significantly re- 
duced but the values between and fc_F become 
possible. 

In fig. 1101 we show the energy transfer distri- 
bution for quasi-elastic events produced with the 
same K2K-like beam. It is clear that it is insensi- 
tive to the effects discussed in the present paper. 
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